The Josephson light-emitting diode 
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We consider an optical quantum dot where an electron level and a hole level are coupled to 
respective superconducting leads. We find that electrons and holes recombine producing photons 
at discrete energies as well as a continuous tail. Further, the spectral lines directly probe the 
induced superconducting correlations on the dot. At energies close to the applied bias voltage eKsd, 
a parameter range exists, where radiation proceeds in pairwise emission of polarization correlated 
photons. At energies close to 2eVBd, emitted photons are associated with Cooper pair transfer and 
are reminiscent of Josephson radiation. We discuss how to probe the coherence of these photons in 
a SQUID geometry via single photon interference. 

PACS numbers: 73.40.-c, 74.45.-|-c, 74.50.-|-r, 78.67.-n 
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Electron-hole recombination in semiconductors accom- 
panied by emission of visible light is a key element of 
many technologies. Semiconducting QDs have been pro- 
posed to enhance these technologies by engineering the 
frequencies of radiation [l[ . In the context of modern re- 
search, they have been considered as a controllable source 
of sins; le i, i, i and entangled two-photon pairs 
QDs allow for integration of photon-based technologies 
and solid state systems where electronic degrees of free- 
dom are used to represent quantum information (e.g. 
electron spins in quantum dots (QDs) [7], charge- [1] and 
flux qubits Q in superconducting (SC) circuits), com- 
bining the advantages of both. For quantum information 
purposes, it is crucial that indistinguishable optical pho- 
tons or pairs of photons can be created on demand. The 
semiconducting QDs provide means to achieve this [10.] . 

SC Josephson junctions can also be a source of coher- 
ent radiation. When the junction is biased with a volt- 
age Vsd^ photons with frequency uo = 2eV^d/^ are emitted 
corresponding to Cooper pair transfers between the two 
SC leads. This radiation is coherent since the Cooper 
pair transfers are coherent owing to macroscopic phase 
coherence of SC condensates involved [Tlj. The frequency 
of Josephson radiation is limited by the SC energy gap 
A ~ 1 meV, huj = 2eVsd < 4A. This is three orders of 
magnitude away from the optical frequency range. 

Many theoretical predictions (e.g. [12]) promote the 
combination of SCs and semiconductors within a single 
nanostructure. This difficult technological problem at- 
tracted attention for a long time [13]. Recent progress 
has been achieved with semiconductor nanowires. SC 
field-effect transistor [14] and Josephson effect [15] in a 
semiconducting QD have been experimentally confirmed. 

In this Letter, we propose and investigate theoretically 
a setup where a superconducting p-n junction enclosing a 
semiconducting QD emits photons in the optical range, 
see Fig. 1. This device is biased by a voltage Vsd which 
is close to the semiconducting band gap. We show that. 




FIG. 1: (Color online) Sketch of a QD in contact to supercon- 
ducting leads with chemical potentials /ie and /jh- We con- 
sider a level in the conduction band (with energy Ee) and a 
level in the valence band (with energy Eh). Each level is only 
coupled to one of the reservoirs as indicated. Photon emission 
processes with energies %oo close to the applied voltage bias 
eVsd and at 2eVsd are indicated. 

owing to SC correlations, the device emits the photons 
in the frequency range eVsd/^ concentrated in several 
discrete spectral lines, the line-width being restricted by 
the emission time only. The acts of photon emission cor- 
relate. In this way, one can arrange emmision of pairs 
of photons of opposite polarization. The device is also 
shown to emit in the frequency range 2eVsd/^' The emit- 
ted light is associated with Cooper pair transfer between 
the SC leads and is therefore coherent. This is in fact 
Josephson radiation at optical frequency. 

Setup details. The semiconducting QD encompasses 
two levels: one for electrons(e), one for holes(h). The 
levels are coupled to corresponding SC leads (source and 
drain), those being characterized by energy gaps Ae,h. 
The levels are aligned to the corresponding chemical po- 
tentials /ie,/i- We count their energies E^^h from these po- 
tentials assuming |^e,^| <^ l^e,/i|- The tunnel coupling 
in the normal state is characterized by the broadening of 
a corresponding level, rt;e,/i, those being proportional to 
squares of the tunneling amplitudes. In the presence of 
superconductivity, we treat the coupling to the SC leads 
in second order perturbation theory [16]. This accounts 
for coherent transfers of electron singlets between the QD 
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and the SC leads, and amounts to an induced pair po- 
tential for the level, with = (1/2) exp[z^e,/i]rt;e,/i 
(assuming Tt-e,h < l^e,/i|, [l2|), and (/)e,/i the phase of 
the corresponding Ae,h. 

The induced pair potential results in formation of four 
discrete low-energy states at each (electron or hole) side 
of the setup. We write the effective low-energy Hamilto- 
nian for electron side, skipping index "e" for A, F, 



D 



■Kc\c\ 



(1) 



where we assume that the charging energy (repulsive on- 
site interaction) U <C |A|. 

By diagonalizing we obtain two degenerate single- 

c||0) with energy 



particle states | t) = c||0) and | [) 



- .t 



E forming a doublet (|0) denotes the empty level), and 
two singlets, those being linear superpositions of |0) and 
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|0). For the ground state singlet, we obtain 
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with energy Eg = E - {E"^ + |A|2)i/2^ {E = E + {lJ/2)). 
The coherence factors ^JJ are |w| = (l/v^)!! ± 
E/{E'^ + |A|2)i/2]i/2^ The excited state singlet reads 
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with energy e^x = E ^ {E'^ + |Ap)^/^. Similary, four 
states are formed on the hole side of the setup. Since 
we are dealing with holes, we define the corresponding 
vacuum \0)h as the level occupied by two electrons |16i] . 
Apart from this difference, the energies and wave func- 
tions of the states are given by above expressions with 
^, A, Ft, /7 = Eh^ A^, Tt-h^Uh' One could easily include 
the interaction energy between electrons and holes in the 
above scheme. We neglect this interaction since we do 
not expect it to change our results qualitatively. 

A SC p-n junction has been discussed in [18], and sup- 
plemented with a QD in [l9], in the context of superra- 
diance which is irrelevant for our proposed effects. 

Emission of '^red^' light. So far, we have not enabled 
charge transfer through the setup. This can only proceed 
by recombination of an electron and a hole at different 
sides of the setup, see Fig. [TJ Such transfer has to dis- 
pose an energy eVsd corresponding the energy differ- 
ence between the electron and hole level, and therefore is 
accompanied by emission of a photon of this energy: Let 
us call it "red" photon. The recombination is described 
by the following Hamiltonian: 

Hint,i = GJ2 (4,-/iiCT +4,+/iTCx) e-^^^-*+H.c. (4) 

The time dependence exp(±iel^d^) accounts for the dif- 
ference between /ig and /i^. Eq. (|4]) is a minimal model for 
the photon-assisted recombination of e-h pairs. We as- 
sume usual selection rules [20] implying that the holes are 
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FIG. 2: (Color online) a) Recombination diagram of the 
biexciton-exciton cascade (EP cycle) with the QD coherently 
coupled to SC leads inducing 4 different singlet states and 
two (degenerate) doublet states for the combined system of 
electrons (e) and holes (h). The cascade produces two "red" 
photons with frequencies ou and ou' (of opposite circular polar- 
ization cr^) on the order eVsd/fi. The cascade can proceed via 
32 different decay paths (illustrated by red arrows) leading up 
to 8 distinct emission peaks as a consequence of induced gap 
Ae,^ in the QD. b) Illustration of "blue" photon emission: 
The biexciton-exciton cascade can also proceed by emission 
of a single coherent photon at the Josephson frequency 2eVsd 
which connects the same initial and final states (e.g. \g)e\g)h) 
differing by the transfer of one Cooper pair. This "blue" 
photon can be "stimulated" by an in-plane dc-electric field 
with Hamiltonian i/int,o- The cascade then involves a "zero- 
frequency photon" (non-radiative decay) as well as a "blue" 
photon and the intermediate state |l)e|l)^ is only occupied 
virtually, whereas the total process conserves energy. 

"heavy" , /i| (/i|) creates a hole with the total angular mo- 
mentum jz = 3/2(— 3/2). Eq. @ then ensures the con- 
servation of total angular momentum: The polarization 
{p = ±) of the photon emitted into a mode q (aj^^) is de- 
termined by the electron and hole spins [21]. An isolated 
QD in the state | t)e| i)h would recombine to |0)e|0)/i 
with the rate Fph cx G^. Since the states of the QD are 
modified by coupling to SC leads [see Eqs. ([2]),([3|)], the 
"red" emission causes transitions between all QD states 
(see also Fig. [2^)). 

An even parity emission cycle (EP) (#e + #/i=even) 
and an odd parity (OP) cycle (#e + #/i=odd) exist. The 
transitions proceed between the discrete states of the QD 
(see Fig. [2^)). They give rise to sharp emission lines with 
frequencies directly related to energy differences between 
the states [22]. The rates incorporate the coherence fac- 
tors, for instance. 



{T^i,/h)\veUh\^; (5) 

{r^^/n)\ueVh\^. (6) 
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EP and OP cycles are connected by transitions of a sec- 
ond type which involve the excitation of a single quasi- 
particle with energy > Ag,/?, in one of the leads [16] and 
therefore change the parity that is conserved in course 
of photon emission. They give rise to a continuous spec- 
trum of the "red" light emitted that is separated from the 
lines by frequency min(Ae, A/j,)/fi. The transition rates 
of the second type are smaller as those of the first type 
by a typical reduction factor |A|/|A| <C 1. 

The emission intensity i{uj) = "^ab p^a^bi^)Pa 
QD can be computed from the probabilities pa to be 
in one of 16 possible QD states \a). They follow from 
the stationary solution of the master equation describ- 
ing the setup dynamics, governed by the rates W^^^ = 
f duj' w^^^^{ijo') [16]. The emission intensity computed is 
shown in Fig. [3] versus photon frequency uo (we assume 
for simplicity that |Ae| = \^h\ and Ue = Uh)- Plot a) 
gives the intensity at the scale \hix) — eVsd\ ^ |A| (for 
the case = = U = 0). Three discrete peaks 
are visible at much smaller scale of the induced gap 
\huo — eVsd\ 1^1- fioo ^ eV^d — |A|, a continuous 
tail of emission starts (enlarged in the inset) reflecting 
quasiparticle creation in the leads. The dashed line is 
the emission spectrum of the same QD without super- 
conductivity. In this case, the spectrum is continuously 
broadened on the scale Ft = 2|A|. The total emission in- 
tensity approximately corresponds to the total intensity 
of the three discrete lines in the SC case. Plot b) illus- 
trates the regime of photon-pair emission. The chosen 
parameters E^ — 1.9, E^ = —1.6 and U = 0.28 (in units 
of I A I) induce a large population of the ground state sin- 
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This has striking consequences for the cascade emis- 
sion process \g)e\g)h ^ |l)e|l)^ ^ \9)e\g)h shown in b) 
(main full lines). From Eqs. (O and (|6]) we deduce, that 

Therefore, this process produces two-photons of oppo- 
site polarization in a pair (i.e. the delay time between 
the emission of the first and second photon is much 
shorter than the emission time of the pair) and with en- 



ergies hco = eVsd ± {Sq 



Ee — Eh). We point out 



that the energies of these correlated photons are differ- 
ent, however, the polarization and energy of the pho- 
tons are uncorrelated. This cascade corresponds to the 
biexciton-exciton decay discussed in [HI in the context of 
polarization-entangled photons. Therefore, potentially 
pairwise entangled photons could be identified effi- 
ciently in the time domain. The dotted lines (OP cycle) 
are energetically distinct from the full lines (EP cycle) as 
a consequence of induced |A| and U. This allows us to 
distinguish emission processes from different cycles. We 
remark that the charge current through the device just 
equals the emission intensity of "red" photons. 

Emission of %lue^^ light at 2eVsd' We now consider 
the emission of a single-photon per Cooper pair transfer 
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FIG. 3: Emission intensity of "red" photons in the energy 
range = \nijO - eVsd\ ^ |A|: We use |A|=0.1|A| and dis- 
crete peaks are broadened with Fph, (Fph/Ft = 0.05 a) and 
0.02 b)). a) Spectrum at resonance Ee — Eh = 0, [/ = 
encompasses several discrete peaks and a small continuous 
tail (inset). The dashed line shows intensity in the case of 
normal leads, b) Regime of pair-emission: Full lines and dot- 
ted lines show the emission spectrum from the EP- and OP 
cycle, resp. Main full lines originate from time- and polariza- 
tion correlated photons emitted from the biexciton-exciton 
cascade with groundstate singlets for electrons and holes. Pa- 
rameters: Ee = 1.9, Eh = -1.6, U = 0.28 (in units of |A|). 

through the QD. Since the Cooper pair charge is 2e, the 
energy associated with its transfer is 2eVsd- If a single 
photon is emitted by this process, it must have a fre- 
quency ~ 2eVsd/h which we call a "blue" photon. Since 
only one electron-hole pair recombines radiatively in the 
emission process, we need a static in-plane electric field 
Eq that annihilates the other pair (see Fig. Wp))- This 
annihilation without emission is described by 
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with (X Eq^x ^ iEo^y. To second order in the total 
interaction Hamiltonian i^int = ^int,i + ^int,0 7 the rate 
to emit a single "blue" photon (with polarization p = ±) 
is = {27r/h)\Al^fS{eb-ea^hu-2eVsd) between 

initial state \a) (with energy Sa) and final state \b) (with 
energy Sb) of the QD. For the case where \a) and \b) 
belong to the singlet subspace [24] , we obtain (in leading 
order in 1/eVsd)^ 



GVo^(6|00){22|a 



2{Ee + Eh) — Eg — £b 



(8) 



We note that the amplitude can also connect different 
initial and final QD states resulting in incoherent pho- 
tons. However, they are emitted at different frequen- 
cies. The light emitted at huo = 2eVsd is always coherent. 




FIG. 4: (Color online) a) Proposed Josephson light-emitting 
diode in a SQUID configuration with two QDs (1,2) enclos- 
ing flux Interference of emitted photons at the Joseph- 
son frequency 2eVsd/T^ is observed via an optical mirror sys- 
tem and photodetector D. b) Intensity Ij of photons at D as 
function of $ (in units of $o = hc/2e) through SQUID for 
Ee — —Eh — |A|, U — in both QDs. We choose equal 
optical path lengths from QDs 1,2 to D. 

In this case, ^^^^ = A^^^ = Al and (a 1 00) (22 1 a) = 
ibexp[i(0e — (l)h)]\ueUhVeVh\' The "blue" emission out of 
the doublet states | t)e| i)h and | |)e| ])h is anomalously 
small with ^^^5 cx (eFgd)"^ and is irrelevant. 

Let us consider two QDs embedded in a SQUID loop 
as shown in Fig. [4^). Coherent emission from either QD 
(1 or 2) into a common photonic mode (and with the 
same polarization) has amplitudes AS^ ^ and A^^ ^, (as- 
suming the QDs are in states |a) and |a'), resp.). The 
total intensity Ij of photons in the common mode is pro- 
portional to Eaa',p PaPa' \A%^^^'-I^' + Al^^,^^^^^ ' \^ , 

where Aj = hcj^eV^^ is the wave length of coher- 
ent light at the Josephson frequency and l\ and I2 are 
the respective path lengths from the QDs to the de- 
tector. The interference contribution is proportional 

to Eaa'.pPaPa'Re[^l„(^l,,)1 with Re[^f )1 CX 

cos[27r((/i — /2)/Aj + <l>/<l>o)]5 where we use that (/)ie — 
<i>\h - (^2e - 02/i) = 27r$/^o, with $ the flux through the 
SQUID and <l>o = /ic/2e the SC flux quantum. 

Fig. IUd) shows the computed emission intensity of 
2el^d photons as a function of flux ^. We find that 
the intensity oscillates with period given by the super- 
conducting flux quantum $0 = hc/2e^ and has a mag- 
nitude of order {V^Y,/n)Q with Q = 2\V^ / (eVsd)^ = 
2\d' Ep\'^ \A\'^ / (eVsd)^ , where d is the optical dipole mo- 
ment of the QD HH. The electric field Eq could 
be created by gates. Typical critical field strengths be- 
fore quenching the photoluminescence of optical QDs are 
on the order of several Volts / //m [2^ . Taking | d 1 1 1 on 
the order of the QD diameter ~ 20 nm and estimating 
I A I < 1 meV (bounded by | A | ) , we arrive at an intensity 
Ij ^ A photons/s assuming eV^d ^ leV and fi/Fph ~ 
0.1 ns at 2eVsd' This intensity is measurable with single- 
photon detectors [27I]. In addition, the Purcell effect in a 
QD-cavity system could enhance Fph substantially [2^. 

In conclusion, we investigated emission from a quan- 
tum dot (QD) embedded in a superconducting (SC) p-n 
junction. The presence of SC leads induces an effective 
pair-potential for electrons (e) and holes (h) on the QD. 
At frequencies uo close to the voltage bias eVgd/^ of the p- 



n junction, a regime exists where radiation is correlated 
in pairs of oppositely polarized photons. At cj = 2eUsd/^, 
emission is associated with Cooper pair transfer and is 
coherent. We proposed an experiment where interference 
of radiation from distant QDs arranged in a SQUID ge- 
ometry can be manipulated by a magnetic flux. This 
provides a fascinating new tool to manipulate coherent 
light at optical frequencies. 
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SUPPLEMENTARY MATERIAL FOR "THE 
JOSEPHSON LIGHT-EMITTING DIODE" 

Spin-degenerate level coupled to a superconductor: 
effective Hamiltonian 

In this section we derive the effective Hamiltonian 
Eq. (1) in the main text for a level (conduction or valence 
band) of the quantum dot (QD) coupled to a supercon- 
ductor (SC). The Hamiltonian for SC and the QD level 
coupled by tunneling is H = Hs -\-Hd ^Ht- The s-wave 
superconductor is described by the BCS Hamiltonian [l*] 

Hs = Y. CkcLck. + E ^kT^-ki + H-c.) , (1) 

kcr k 

with = Sk — l-i the single-particle energies in SC 
counted from the chemical potential ji. This Hamil- 
tonian is diagonalized by the canonical transformation 

CkT ^fc7kT + ^/eTlki ^-ki ^fc7-ki - ^/e7kT 

reads Hs = Eka^^TkaTka with Ef, = y^^l + |Ap, and 
A the superconducting pair-potential. The isolated QD is 
represented by Hd = ^Eo- ^t^(T~^Un^ni, {a =t, |) with 
U a possible repulsive on-site interaction and E the spin- 
degenerate energy level (counted from fi). The tunnel- 
ing Hamiltonian has the form Ht = Ekcr ^kcJ-Cko- + H.c. 
The first step is to integrate out the SC by deriving an 
effective QD Hamiltonian in the subspace of the BCS 
groundstate taking into account the tunneling between 
the QD and SC. Defining P as the projection operator 
for states of the total system with no excitations in SC, 
i.e. jkaP'^ = for any state ?/^, the effective Hamiltonian 
isi 

Hd{s) = PHP + PH \^QHP, (2) 

with Q = 1 - P. The first term on the RHS of Eq. ^ 
is replaced by Hd since the tunneling Ht cannot act to 
first order in the subspace with projector P. To second 
order in Ht^ the resonant transport of electron singlets 
between the QD and the SC is possible and described by 
the second term on the RHS of Eq. ([2]). To leading order 
in Ht, this gives PH{e-QH)-^QHP = PHT[e-{Hs^ 
Hd)]~^HtP. The virtual energy cost e - Q{Hd + Hs) 
(created by hopping of a single electron from (to) the QD 
to (from) the SC) with a quasiparticle of energy Ek in 
the SC is approximated by —Ek since we assume that 
A is the largest energy scale, i.e. |A| \E\,U,tk. By 
tunneling of another electron (with opposite spin) from 
the QD to the SC (or vice versa), the excitation in SC 
can be removed and a Cooper pair is added (or removed) 
to (from) the condensate. These processes lead to the 



following contribution 

PHT[e-{Hs^HD)]-^HTP 

--El^k|'^(^i^T-CTCi)+H.c., (3) 
k ^ 

where we used that = t\^. Since we are dealing with 
two SCs (electron-side and hole-side of the setup), it is 
important to keep track of the SC condensate phase <p 
which is related to the phases of the SC coherence factors 
\U' <^k = -\vk\uk\ex.p{-i(j)), where \uk\ = (1/V2)(1 + 
a/^fc)'/', \vk\ = (l/y2)(l-a/^fc)^/^ By replacing the 
momentum sum in Eq. (|3]) by an integral over energy, we 
obtain the following effective QD Hamiltonian 

Hd=eY^ c\c^ + Acjcj + A*qc| + Un^n^ , (4) 

with A = (1/2) exp(z0)rt, where the level broadening 
Ft = 2i\vs\t\^ ' Here, vs is the normal-state DOS per spin 
at the Fermi level \i in SC. Since the main part of the in- 
tegral in Eq. (|3|) comes from energies \^k\ ^ |A| <C (^f 
the Fermi energy in the leads), the exact /c-dependence 
of tk can be neglected. 

Energy levels of the QD coupled to SC leads 

The diagonalization of Hb leads to four states for elec- 
trons and holes, see Fig. 1. For the electron side of the 
setup (with bare level energy E^.^ induced gap Ag and 
on-site repulsion Uq) there is one doublet state 

|T)e = c||0)e (5) 

and 

U)e = c||0)e, (6) 

with energy £^e, and two singlets (being a superpositions 
of zero and two electrons) 

|5)e = -e-'*«|Ue||0)e + |Ve||2)e, (7) 

with 

e% = K- \jEl^\K\\ (8) 

where = -^e + f/e/2, |2)e = c|c||0)e and |0)e de- 
notes the empty level. We have introduced the coher- 
ence factors |We| = {\l + E,l{El + | Ae|2)i/2]i/2 
and \ve\ = (1/V2)[1 - ^e/(4' + |Ae|2)i/2]i/2. The ex- 
cited state involving the superconductor is 

|e)e = e-'*«|«e||0)e + |We||2)e, (9) 

with 

el, = Ee + ^El + \K?. (10) 



6 



a) 



0(2) ^ex 
— — — — —^^^^^ 


2{qV 




e-W 0(2) 


-4 -2 2 4 



&) 

0.6 

^ 0.2 

0.0 



li 



WW 



IJ 



M?/|A| 



FIG. 1: a) Level structure of the QD coupled to SC leads: 
The coupling to SC leads induces four states and three dis- 
tinct energies for electrons (e) and for holes (h). A doublet 



with energy Ee(h) and two singlets with energies Sg^^'' and 
Sei^'* being superpositions of zero and two electrons (holes). 
The mixing character can be tuned by the QD level en- 
ergy Ee(h)- Away from resonance, |^e(^)| ^ the states 
become pure number states (shown with numbers for elec- 
trons(holes), energies given by dashed lines), where A is the 
induced SC gap in the QD. The level on-site energy is cho- 
sen as [/ = 0.1|A|. b) Emission intensity of "red" photons 
at hcj ^ eVsd- In the generic case shown (Ee = —0.81, 
Eh = -0.53 and U = 0.47 (in units of |A|)) 12 peaks are 
visible (assuming |Ae| — |A^| = |A| and Ue = Uh = U). 
Peaks are broadened with Fph (Fph/Ft = 0.02). 




FIG. 2: (Color online) Setup dynamics and "red" photon 
emission: QD number states and transitions between them in- 
volving emission of "red" photons: Numbers within circles de- 
note occupation numbers for electrons (left) and holes (right). 
Full lines (red) depict transitions (direction given by the ar- 
row) where the total (e+h) number of particles on the levels 
change by two. Dashed black lines denote coherence between 
states mixed by SCs. Two such cycles exist: The upper one 
involves only states with an even number of total particles 
(even parity (EP) cycle) whereas the lower cycle has an odd 
number of total particles (odd parity (OP) cycle). The two 
cycles are connected by processes (dotted lines) which create 
a quasiparticle excitation in one of the leads (brown for n- 
side, green for p-side) in combination of emission of a "red" 
photon. 



For the hole-side of the setup the same four levels result 
(with Ee Eh, Ae Ah, and Ue Uh, ^ M- We 
then transform to the hole-picture for the valence band, 
by defining \0)h = |2)e and Ca = h^-a- Explicitly, the 
four levels on the hole-side are 

\^)h = h\\0)n (11) 

and 

\i)h = h\\0)h, (12) 
\g)h = -e-"^'\uh\\2)h + \vh\\0)h (13) 

and 

|e)/, = e-**^|«^||2)^ + |«„||0)^, (14) 
with 12),. = h\hl \0)h. 



"Red" photon emission rates and master equation 

In this section we derive the rates for spontaneous 
emission of photons due to electron-hole recombination 
which is also responsible for charge transport through 
the QD. The interaction Hamiltonian with the elm- field 
takes on the form (Eq. (4) of main text) 

Hint,i =GY, +4,+ ^TCi) e-*'^^^' + H.c. 

q 

(15) 

The circular polarization (jp = ±) of the photon emitted 
into a mode q {a^^^) is determined by the electron and 
hole spins [3], and the Hamiltonian for photons is i^ph = 

We treat the interaction Hamiltonian with the elm- 
field i^int,i as a perturbation and assume Fph <C |Ae./i|, 
where Fph = 27ri^ph|Gp with i^ph the photon DOS per po- 
larization direction, assumed to be independent of energy 
and polarization. We use a stationary master equation 
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approach to calculate the occupation probability pa for 
each of the 16 possible states \a) of the combined system 
of electrons and holes. These states \a) are: \g)e\g)hi 

|e)e|e)h, \g)e\e)h, \e)e\g)h, \g)e\ T)/i, \g)e\ i)h, |e)e| T)h, 
\e)e\ I ])e\g)h^ I [)e\g)h. \ ])e\e)h, \ [)e\e)h, \ ]) e\ 

U)e| T)/., I Del T)/., IDeli)/.. 

The states are connected by rates of the form W^^, 
where \a) is the initial state of the QD (energy £a) and \h) 
is the final state of the QD (energy £5) via the emission of 
a photon with energy Tiuj and polarization p = ±. They 
are given by the usual form Q 

= y ^ l(&;9,p|^fint,il«;0)|' 5{sa - Sb - hcb,), 

(16) 

with uj = uj — eVsd/^- The dynamics of the system is 
illustrated by the diagram in Fig. 2. Two emission cy- 



cles exist: A cycle where #e + is even (upper cycle) 
which we refer to as the even parity (EP) cycle and a 
cycle where #e + #/i is odd (lower cycle) which we refer 
to as the odd parity (OP) cycle. Full red lines connect 
states within the same cycle (in the direction of arrows) 
via emission of a photon with energy eV^d- The two 
cycles are connected by rates of a second type (illustrated 
in Fig. 3 for a specific example). Since the parity cannot 
change in the course of photon emission only, these cycle 
connecting processes (depicted by dotted lines in Fig. 2) 
create an excitation in one of the leads via single-particle 
tunneling. Since eVsd ^ ^e,h such processes are pos- 
sible in combination with photon emission and will be 
discussed in more detail below. 

Formally, the dynamics is governed by the following 
master equation (we abbreviate the QD states \a)e\P)h 
as aP) 



W+ +W|" +Wt +W~ +W~ 

U,gg ^ Ihgg ^ Tg,gg lg,gg ghgg ^ gl,gg 



Pgg, (17) 



Pee = W^e.n PU + ^et.iT ^iT + Ke,U PU + W^etie Pie + W+^e]' Pel + Ke,el Pel 

- Wfuee + W^fT.ee + K,ee + W[e,ee + We'^ee + W^ee] Pee, (18) 



Pge = Wge.n PH + ^/e.iT ''iT + W^ge.Te Pie + W+^le Pie + W+^g^ Pg^ + W^^^^^ Pg^ 



w+ +wr, +W+ +wr +w~, 

ll,ge ^ II, ge ^ le,ge ^ " le,ge ^ " gl ,ge ^ " gl,ge 



Pge, (19) 



Peg = W^eg.n Pi I + ^^g,ll PH + ^egAg Pig + W^g,l9 Pl9 + ^tg,e\ Pel + W^eg.ei Pel 



Wll,eg + W-,^eg + W^g,eg + W-^,eg + W-^eg + W:,^eg Peg, (20) 



Pgl = ^aT.To Pig + ^,T,Te PU + ^ „^ ,„„ Pgg + W^ot.ge Pge + W^,"t, I T Pll + W^„T,TT ^TT 



Tg.sT ' le,g] ' ' gg,gl ' ge,g] ' ' iT.sT ' TT.sT '^fl' ^ ' 



P9i = + Kide Pie + ^ah.. P99 + Kl.c Poe + ^+i,U PU + ^ahTl PU 



gi,ge ' 



gidi 



W[g,gi + W^x"e..i + W-,g, + W^,-,,, + W^A.i + ^U.i 



pgl, (22) 



Pel = W"eT,Ts P^9 + W^eT.Te Pie + ^^eT.ee Pee + Pll + ^etiT Pll 

- K,e, + ^eteT + K,e, + K,el + K ,el + Pel, (23) 



Pel = Wtue Pie + W+^lg Pig + W+^ee Pee + W^egPeg + W^etii PH + ^ei.Ti Pll 



Pel, (24) 
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hg = Wt,. PeT + W+ pg^ + W+ pgg + W+ Peg + W+ ptT + W... pn 



^."t,T. + ^el,T. + W-,,g + W--^,T. + ^U,T. + ^n.T. 



Pla, (25) 



Pi9 = ^l9M Pel + W^g,glP9l + ^l9,99p99 + ^I9,e9pe9 + W^g,UPll + ^itiT^'iT 



Kug + K^g + K^g + K^ig + W^-Ug + Kug Pl9, (26) 



PU = W^tt.sT P91 + ^T^e.eT Pel + W^e,ge P9e + W^e,ee Pee + W^t+TT ^TT + ^T^.Ti Ptl 



PU, (27) 



Pie = W^e,9l P9l + W^e,el Pel + W^e,9e P9e + W^e,ee Pee + M^ie.ii PU + W^e,n Pit 



Kue + Kue + WeXie + K,le + ^^^.e + W^Ue 



Pie, (28) 



Pn = W^l,gg P99 + W^l,ge Pge + W"u,es Pe9 + W^U,ee Pee + W^U,Te P^e + ^fugl Psl + ^U,ei Pel + ^U,Ts ^Ts 



PTi, (29) 



Pll = W^T,ie Pie + W^T,ig Pig + W^iT.gT P9I + ^iT.eT Pel + ^Ihgg Pgg + ^Ihge Pge + W^^,eg Peg + W^^,ee Pee 

- n,iT + K,ii + Ki,ii + K,ii +K,ii + K,ii + Ke,u + pii' (30) 



Pll = ^TT.TgPTs + ^TT.TePTe + W^y,g^Pgl + W^^,e-[ Pel 



^tXtT + ^T9,TT + ^9T,TT + ^eT.TT 



Plh (31) 



Pll = ^lUgPlS + ^lUePle + W^l,glP9l + ^IMPel 



w[e,u + ^r.,u + Kui + Kui 



Pll- (32) 



r 



The transition rates within the same cycle have the fol- 
lowing form, e.g. 



^ll,gg = {^pi^/^)\VeUh\\ 



(33) 



which emits a (7+ -photon at energy htu = Sg -\- Sg — Eq 
Eh. 

An example for the OP cycle is 



(34) 



which emits a (T+-photon at energy Tiw = s^^ —Sg + Ee — 
Eh- 

The rates that connect the two cycles involve the tran- 
sition operator V{sa - Ho)-'^V with Hq = Hf^ -\- -\- 
Hs + i^ph, V = i^int,! ^ and Sa is the energy 

of the QD before the transition (i.e. \a) is an eigen- 
state of Ho with no quasiparticle in SC leads and no 
photons present). These rates are different since they in- 
volve the tunneling of an electron (hole) into/from the 



SC reservoirs (creating a quasiparticle with energy of at 
least |Ae,/i|), in combination with emission of a photon 
such that the total energy is conserved in the final state. 
These processes have the following rates, e.g. 



w 



S {nCj -s^^Ee^ E{Ck)) , (35) 



which creates a cr+-photon at energy hoj = Sg — Ee — 
E{^k) and a quasiparticle in the SC reservoir (hole- 
side) with spin down, momentum k and energy E{£^k) = 
V^K|A^|2, i.e. the state 7LilO)BCS, see Fig. 3. We 
can integrate over the quasiparticle state in the lead to 
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get the rate for photon emission at frequency uj 



w 



^9,99 



(u) Tph 



Tiuo) 

- Ee — hu 



|A.|) 



(36) 



Note that the spectrum of "red" photons contains a con- 
tinuous tail (see Fig. 3a) of main text) due to these pro- 
cesses. For the master equation, we need the total emis- 



sion rate W^g 



^9,99 



J dujw^g 99^^^^ with the result 



(Fph/^) \Ve 



2Ah 



(37) 



We remark that these cycle connecting rates also al- 
low the population of triplet QD states, | t)e| T)/i and 
I i)e| i)h and within our model are also responsible for 
the decay of triplet QD states that cannot proceed by di- 
rect recombination owing to selection rules (see Eq. (15)). 
Similar results as in Eqs. (|33]), (|34]), (|36]) and ^ hold 
for all processes that are included in the master equation 
[Eqs. (17)-(32)]. In the cycle connecting processes, we 
only include terms to leading order in powers of l/Ae^h- 



Spontaneous emission of "red" photons 

We solve the master equation (Eqs. (17)- (32)) in the 
stationary limit pa = and calculate the emission inten- 
sity according to 



ab,p 



(38) 



which leads to the plots of Fig. 3 in the main text, and 
to Fig. lb) in this supplementary. 

We now discuss the case of normal leads. When the 
leads are normal conducting, single particle transfer be- 
tween the leads and QD is possible and the QD levels 
acquire the usual broadening (Ftje for the electron level 
and Ft;^ for the hole level). Since we are working in the 
parameter regime Tt-e^^t-.h ^ Fph, the level broadening 
due to the electron hole recombination can be neglected 
here. This then leads to the following emission intensity 
(for = 0) at frequency uo 

+ CO +CO 

iN(w) = 2rph [ de [ de'fe{e)iye{e)[l- fh{e')] 



-CO — OO 



X iyh{e')S{e-e' -hco), (39) 

where fe,h{^) = (1 + exp[/3(£: — iie,h)])~^ are the usual 
lead Fermi functions with (5 = {kBT)~^ and Ve,h = 



Q:(20) + /3 (22) 




(2J)U @+hv 




FIG. 3: Cycle connecting processes: Sketch of the cycle con- 
necting process (proceeds in direction of arrows) that switches 
the QD state from the EP cycle with electron and hole singlets 
to I ])e\g)h or I ])e\e)h from the OP cycle (see also Fig. [2]) 
and which creates a quasiparticle (7^^^^) in the SC lead (hole- 
side) with momentum k and spin [ thereby emitting a "red" 
photon with frequency v and circular polarization +. Only 
the initial state component with two electrons will take part 
in the process, whereas both components of the hole singlet 
(|0)^ or 1 2)^) can participate in the process. The coherence 
factors a^[3^a\[3' depend on specific QD states (initial and 
final states) involved. This diagram is therefore also relevant 
for the transition \g)e\g)h \ '\)e\g)h discussed in the text. 



(l/^)(rt;e,/,/2)/[(£ - Ee^H + Pe^? + (rt;e,/^/2)2] are the 
DOS for the broadened QD levels. The factor 2 in 
Eq. (|39|) accounts for the two polarizations of the emit- 
ted light. In Fig. 3a) of the main text we show (dashed 
curve) i]<i{uo) at zero temperature (T = 0) and at reso- 
nance Ee = Eh = with the result (assuming Ft;e = 
Et;^) iNioo) = (Fph/|A|)(2/^2)^(^~/|^|) ^^^^^ ^ 



Q{-z) /q ^ dx[{x^^l){{x^z)^^l)]-^. We note that the 
three discrete peaks in Fig. 3a) of the main text (SC case) 
approximately corresponds to the integrated emission in- 
tensity in the normal case and the leading contributions 
for large negative uj (continuous contribution in SC case) 
fall off in both cases like l/Co'^. 



Spontaneous emission at the Josephson frequency 

2eVsd/n 

Here, we describe the process of photon emission at 
the Josephson frequency 2eV^d/^- If only one photon 
per Cooper pair transfer from the n side to the p-side is 
emitted it must have an energy ~ 2eVsd which we call a 
"blue" photon. Since the Cooper pair has charge 2e one 
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electron-hole pair has to recombine without the emission 
of a photon which becomes possible in the presence of 
an externally applied dc electric field Eq. We now derive 
the form of the relevant Hamiltonian Eq. (7) of the main 
text. 

Besides the Hamiltonian i^int,i that emits or absorbs 
photons, there is an additional part of the Hamiltonian 
related to a dc-electric field Eq 



^int,o = -er- Eo, 



(40) 



where we assume that the field is homogeneous. Since 
the field is static it cannot provide photons. In second 
quantization, i^int,o reads 



(41) 



and 



smt 



(44) 



To calculate Vq ^ ^ we make use of the fact that the 
envelope part of the wave functions vary slowly on the 
scale of the lattice and write r = Vi -\- Ri with Ri the 
Bravais lattice vector of the i-th unit cell p\. Using the 
orthogonality of the periodic parts of the wave functions 
from the conduction band and valence band, we obtain 
V^'"" = V^5a',a with = da-Eo. Explicitly, the inter- 
band dipole moment of the QD is da = d{aex — iCy) / a/2, 
where 



d=^Y.^l{R^)^c{R^) I drr'R:{r)R,{r), (45) 



with = {0\ba'{-erEo)cl\0). Here, and de- 

note annihilation operators for electrons with spin a for 
the level in the conduction band and valence band (their 
energies are again counted from respective chemical po- 
tentials jiie,h)^ respectively. 

To calculate the matrix element Vq ^, it is crucial to 
know the orbital angular momentum of states that are 
connected by the operator Eq. (|4Q|) . For QDs, usually the 
conduction band ground state level is an s-state (/ = 0) 
and the valence band state is a heavy hole p-state (/ = 1). 
Both states are two-fold degenerate with opposite total 
angular momentum jz in z-direction. For the conduc- 
tion band these are |l/2, +1/2) and |l/2, —1/2), and for 
the valence band they are |3/2,+3/2) and |3/2,— 3/2). 
Note that valence band states with either total angu- 
lar momentum 3/2 but projection ±1/2 (light-holes) or 
the split-off band with total angular momentum 1/2 are 
lower energy states and are therefore occupied and far 
away from resonance with the leads. 

The wave functions for the s-state in the conduction 
and for the heavy-hole (HH) p-state in the valence band 
are written in the envelope approximation as (r|cj.|0) o:^ 
(l)c{r)uca{r)\cr) and {r\bl\0) c:^ (l)v{r)uuua{r)\cr) , respec- 
tively. Here ^c,v are the envelope functions for the elec- 
tron level and hole level, respectively, and Uca{r) and 
'^HHcr(^) are the k = "Bloch" -parts of the wave func- 
tions which have the periodicity of the lattice. They 
refiect the symmetry of the band and have the form 
UcAr) = Rc{r)Yo°{e,^)\a), unH^{r) = R,{r)Y+\0,^)\ T 
) and iiHHi(^) = Rvi'^)^! "'^(^7 0)1 i)- The spherical har- 



monics are 



e** sin( 



(42) 



(43) 



with i?w the radius of the Wigner-Seitz cell (assumed 
to be sperical for definiteness). The amplitude of the 
dipole moment Eq. (|45|) depends on the specific form 
and material of the QD. To transform to hole opera- 
tors we replace ba- by h^_^ in Eq. (|4T]) which leads to 
Eq. (7) of the main text. Having clarified the process of 
electron-hole recombination without photon emission we 
now can calculate the emission rates for emission of a sin- 
gle "blue" photon, treating the total interaction Hamil- 
tonian i^int = i^int,o + ^int,i as a pcrturbatiou. We first 
discuss the emission of a "blue" photon out of the EP 
cycle (see Fig. 2 upper cycle). Starting in the singlet 
subspace of electrons and holes, the emission process can 
be considered as a second-order "emission" process where 
the subspace of the doublet states {| t)e| i)/i,| i)e| T)/i} 
can only be occupied virtually, and the real transition 
therefore connects the QD singlet subspace with itself 
and produces a photon of hco 2el^d since a charge 2e 
(in form of a Cooper pair) has been transferred by this 
process from the n-side to the p-side of the setup. Note 
that one of the electron-hole pairs present in the initial 
state is annihilated by i^int,o and the other pair by the 
spontaneous emission of a "blue" photon. The annihila- 
tion due to the dc-field can therefore also be interpreted 
as a stimulated emission of a zero frequency photon. 

To second-order in time-dependent perturbation the- 
ory in i^int, the emission rate of a "blue" photon with 
polarization p has the usual form Q 



27r 

T 



^ (/|i^int|^)(^|^intK) 
m 



- Em^ eVsd 
X 8{Ef - E, 



2eFsd), (46) 



where |z), |m), and |/) are the initial, intermediate and 
final states of the QD (energies counted from the respec- 
tive chemical potentials) and elm-environment (photons) 
in the absence of the electric field. Note that the ex- 
plicit time-dependence of i^int is accounted for by the 
el^d-terms in the rate Eq. (|^ . 
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FIG. 4: (Color online) Emission intensity of proposed Joseph- 
son light-emitting diode in a SQUID configuration with two 
QDs (1,2) [see Fig. 4a) of main text] enclosing flux a) 
Generic (asymmetric) situation with Ee = 0.2, Eh = 0.18, 
U = 0.25 in QDl; Ee = 0.1, Eh = -0.3, U = 0.15 in QD2. 
b) Dependence of interference part of Ij (Pf^) at $ = as 
function of Ee and Eh in QDl fixing Ee — 1, Eh = — 1 for 
QD2 (U — in both QDs). c) Regime where QDs are mainly 
in the state of ground state singlets \g)e\g)h with Ee = 2.8, 
Eh = -3.1, Ue = 0.23 in QDl, Ee = 2.9, Eh = -3.2, Uh = 0.6 
in QD2, where almost destructive interference can be reached. 
For plots we choose |Ai;e,h,| = |A2;e,^| = equal optical 
path lengths from QDs 1,2 to D, energies in units of |A| and 
define Q = 2|yo^A|V(eVsd)". 



We are interested in the emission of coherent photons, 
i.e. emission processes that leave the QD state unchanged 
after the photon is released. Such photons are emitted 
at the Josephson frequency 2eVsd/h. Therefore the rate 



Eq. (|46|) can be written as W 



f4 



(rphM)I^SP, 

where a denotes the initial (and final) state \a) of the 
QD. We derive the following amplitudes: 



4) 



sd) 



(47) 



2(J?e+i?/»-gL-geJ 



(eT4d)2 



(eT4d)2 



(48) 



(49) 



(50) 



We note that emission processes out of the doublet 
states |cr)e| — cr)h turn out to be of higher order (cx 
l/(eVsd)^) and will be neglected. We defined = 

V^G\UeVeUhVh \ exp[z(^e - ^/i)]- 

Within the OP cycle, similar amplitudes exist: 



Al„ = AP 



(ev; 



sdj 



and 



AP 



-AP 



Ml 



(eysd)2 ■ 



(51) 



(52) 



(53) 



(54) 



As can be seen from Fig. 2, the "blue" photon emission 
out of the OP cycle involves virtual QD states |m) that 
are a product of a doublet state and a singlet state. 

To calculate the emission intensity at the Joseph- 
son frequency /j, we use the same master equation 
[Eqs. (17)- (32)] since the occupation probabilities of 
the QD states are determined by the faster "red" 
photon emission (see discussion below), i.e. Ij = 
v^^^ci^'^'^ ^^^^ coherence of such pho- 

tons we suggest an interference experiment of photons 
emitted from either of two QDs arranged in a supercon- 
ducting quantum interference device (SQUID), see Fig. 4 
a) of the main text. The emission intensity at the de- 

tector I) is /j = (^ph/^)Eaa^pPaPa'|^?,ae^'"^^/^^ + 

^f^,e^2^^2/Aj|2^ where Aj = hcj^eV^d is the wave length 
of coherent light at the Josephson frequency and l\ and I2 
are the respective path lengths from the QDs to the de- 
tector. The interference contribution Ff^ is proportional 

to Y.aa.,,PaPa'^^\AlMla'r\ with Re[^t, (^^,„, )1 CX 

cos[27r((/i — /2)/Aj + ^/^o)], where we use that — 
<i>\h - (^2e - ^2/1) = 27r^/<l>o, with $ the flux through the 
SQUID and <l>o = /ic/2e the SC flux quantum. In Fig. 4b) 
of the main text we show the emission intensity in a 
regime where we observe a maximal interference contri- 
bution. Note that exactly at resonance = Eh = U = 
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in both QDs, the interference contribution vanishes due 
to different signs in amphtudes [see Eqs. (47)- (54)]. The 
interference contribution can be on the same order as the 
total emission intensity Ij for a quite general parameter 
set [see Fig. H^)] and has the order of magnitude 

ij ~ 2(rphA)|d • £;o|2|A|V(ev;dr- (55) 

Therefore, the "blue" photon emission intensity is ap- 
proximately by a factor \d- Eq\^\K\^ / {eV^dY smaller than 
the "red" photon emission. We note that the bias volt- 
age eVsd over both QDs in the SQUID is necessarily the 
same. This leads to certain constraints regarding the sim- 
ilarities of the two QDs. Gate voltages, however, could 
be used to tune the QD levels in the conduction and va- 
lence band into the close proximity to the SC reservoirs. 
Fig. [4)3) shows the sensitivity of the coherent contribu- 
tion (/}"*) to the change of electron and hole energies of 
one of the QDs (leaving the parameters of the other QD 
fixed). The plot shows that the interference contribution 
Pf^ changes on the scale of |A| and therefore is not very 
sensitive to spectroscopic differences of the two QDs. In 
particular, the two QDs need not be identical within the 
linewidth Fph <C |A|. This linewidth (or broadening) 
of the 2eVsd emission is determined by the much faster 



"red" photon emission which switches between different 
QD states. Fig. [It) shows the regime where the QDs 
are mostly in a single state {\g)e\g)h for both QDs) and 
almost destructive interference can be reached. Finally 
we note that the present phenomena are inherently of SC 
(Josephson) origin, since /j = if \ueUhVeVh\ = which 
is the case if Ag = or A^ = 0. 
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